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o : 

, We study the 2D Ising model on a square lattice with additional non-equal diagonal next-nearest 

, neighbor interactions. The cases of classical and quantum (transverse) models are considered. Pos- 

sible phases and their locations in the space of three Ising couplings are analyzed. In particular, 
^ I incommensurate phases occurring only at non-equal diagonal couplings, are predicted. We also ana- 

lyze a spin-pseudospin model comprised of the quantum Ising model coupled to spin chains in a 
particular region of interactions, corresponding to the Ising sector's super-antiferromagnetic (SAF) 
CO ' ground state. The spin-SAF transition in the coupled Ising-XY model into a phase with co-existent 

SAF Ising (pseudospin) long-range order and a spin gap is considered. Along with destruction of 
the quantum critical point of the Ising sector, the phase digram of the Ising-XY model can also 
demonstrate a re-entrance of the spin-SAF phase. A detailed study of the latter is presented. The 
^ , mechanism of the re-entrance, due to interplay of interactions in the coupled model, and the con- 

. ditions of its appearance are established. Applications of the spin-SAF theory for the transition in 

the quarter- filled ladder compound NaV2 05 are discussed. 

PACS numbers: 71.10.Fd, 71.10.Hf, 75.30.Et, 64.60.-i 
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S ■ The role of competing interactions in ordering is a fascinating problem of condensed matter physics. One of the most 
' canonical examples of such systems are frustrated Ising models which demonstrate a plethora of critical properties, far 
(~| from being exhaustively studied. (For a review see The frustrations^ can be either geometrical, like, e.g., in the 
O Ising model on a triangular lattice, or they can can be brought about by the next-nearest neighbor (nnn) interactions. 
, O, , Competing interactions (frustrations) can, e.g., result in new phases, change the Ising universality class, or even 
destroy the order at all. Another interesting aspect of the criticality in frustrated Ising models is an appearance of 
P\J , Quantum Critical Point(s) (QCP) at special frustration points of model's high degeneracy, and related quantum phase 
^ transitions 

\ Inclusion of a transverse field (f2) brings an extra scale into the game, giving raise to a new and complicated critical 
. behavior. The Ising models with f2 = and Q ^ are also often called classical and quantum, respectively. For a 
' review on the Ising Models in Transverse Field (IMTF) see "S]. Most studies of the frustrated quantum Ising models 
are restricted to their ground states properties, when mapping of the d-dimensional quantum model at T = onto its 
(rf + l)-dimensional classical counterpart helps to analyze the ground state phase diagram of the former. For the nnn 
2D models we are interested in, there has been a considerable effort on the quantum ANNNI model, reviewed in 
The studies of some other 2D frustrated transverse Ising models have appeared only recently 0, . 

Our interest in the subject comes from the earlier work on a quantum Ising model coupled to the spin chains |^. 
H This kind of coupled so-called spin-pseudospin (or spin-orbital) models appear in context of the phase transition in 
NaV205, which has inspired a great experimental and theoretical effort in recent years. (For a review see 0.) Going 
deeper into analysis, we came to realize that the Ising sector of the problem is in fact the 2D transverse nn and nnn 
Ising model on a square lattice. It turns out that even its classical counterpart (fi — 0) was studied only for the 
case of equal nnn couplings Ji — J2 Q ■ The elementary plaquette of this lattice with the notations for couplings is 
shown in Fig. ^ To the best of our knowledge, this 2D nn and nnn Ising model in transverse field is a complete terra 
incognita even at Ji = J2- 

■ So, as the first step, we find the ground state phase diagram of the classical (f2 = 0) Ising model at arbitrary 
5—1 ' couplings Jn, Ji, J2. Along with the three ordered phases found earlier by Fan and Wu I3| for the case Ji = J2 
-ferromagnetic (FM), antiferromagnetic (AF), and super-antiferromagnetic (SAF)- the model has a fourth phase 
which can occur if sign(Ji/J2) = —1. We call it super-ferro-antiferromagnetic (SFAF) j^. From mean- field- type 
arguments we predict also the existence of an incommensurate (IC) phase at T > in this model. The 2D IC phase 
is also called floating ^3 • Similar phase is known for the well-studied 2D ANNNI model P, 0, ^ . We present a 



^ We use the term "frustration" in a broad sense meaning only that there is no spin arrangement on an elementary plaquette which 
can satisfy all bonds. 
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qualitative temperature phase diagram for the regions of the coupHng space where the IC phase is located. Note that 
the three phases - SAF, SFAF, IC - can occur only in the presence of competing interactions in the Ising model, and 
the latter two occur only if Ji ^ Ji. 

This analysis of the classical Ising model lays the grounds for venturing into its study in presence of a transverse 
field. The role of transverse field is subtle. A more straightforward aspect is that its increase above certain critical 
value can eventually destroy the ordered state of the classical model, and in the ground state the transverse field 
results in appearance of a QCP. This is similar to the well-understood quantum nn Ising model. Another particularly 
interesting aspect in the role of transverse field is that it can lift degeneracy of the ground state and stabilize new 
phases at finite temperature in a highly frustrated model, like, e.g., the antiferromagnetic isotropic triangular Ising 
model 13,0], which is disordered at any T > when = 0. 

The behavior of the systems with infinitely degenerate ground states (with or without a finite ground-state entropy 
per spin) can be quite complicated in the presence of transverse field. It lies beyond the scope of the present work, 
and definitely cannot be understood from the mean-field analysis we apply in this study. For the nn and nnn Ising 
model we only identify the lines (planes) in the space of couplings (Jq, Ji, J2) where the model is highly degenerate, 
and in their neighborhood we expect some new exotic phases generated by f2 7^ to appear. 

From mapping of the nn and nnn IMTF at T = onto its classical 3D counterpart, we qualitatively predict the 
(mean- field) ground-state phase boundaries of the quantum model in the coupling space (Jn, Ji^J-i). In particular, it 
follows from our analysis that in the presence of transverse field the IC ground-state phase can penetrate into some 
parts of the FM, AF, SAF regions of the classical model (fi = 0). 

Finally, we consider the coupled spin-pseudospin model. It is proposed to analyze the transition in NaV205. This 
material provides a unique example of a correlated electron system, where the interplay of charge and spin degrees 
of freedom results in a phase transition into a phase with coexistent spin gap and charge order. NaV2 05 is the only 
known so far quarter-filled ladder compound. Each individual rung of a ladder is occupied by single electron which 
is equally distributed between its left/right sites in the disordered phase. At = 34K this compound undergoes a 
phase transition when a spin gap opens, accompanied by charge ordering 0. 

The problem of the electrons in NaV2 05 localized on the rungs of the 2D array of ladders is mapped onto the coupled 
spin-pseudospin model on the effective square lattice. The Ising sector of this model is given by the Hamiltonian of the 
nn and nnn IMTF, and the Ising variables (called pseudospins for this case) represent physically the charge degrees of 
freedom. We model the spin sector by the array of the XY spin chains. The 2D long-range charge order in NaV205 
is identified as the SAF phase of the Ising model, and we restrict our analysis to the SAF region of the couplings 
(Jn, Ji, J2). The coupled model is handled by combining the mean- field treatment of its Ising sector with the use 
of exact results available for the XY spin chains. Since the SAF state is only four-fold degenerate, the mean-field 
predictions for the Ising sector of the coupled model are expected to be at least qualitatively correct. 

The mean-field equations for the coupled model are, with some minor modifications, the same as we have obtained 
earlier |^. A striking feature of the coupled model is that it always orders from the charge-disordered spin-gapless 
state into the phase of co-existent SAF charge order and spin gap. We call this the spin-SAF transition. By always 
we mean that the critical temperature of the spin-SAF transition is non-zero for all Ising couplings within the whole 
considered SAF region. In other words, the QCP of the IMTF is destroyed, and this is due to the spin-charge 
(-pseudospin) coupling. This property of the spin-SAF transition and the parameters of the spin-SAF phase were 
studied earlier so in this work we only reinstate some points and stress the distinctions pertinent to the present 
model. 

The other remarkable feature of the coupled model's phase diagram is re-entrance, which was not well understood 
in our earlier work Q . Now we carry out an analytical study of the re-entrance and establish the conditions when it 
can occur. This analysis allows us to understand the detailed mechanism of this interesting phenomenon generated 
by competing interactions. 

The rest of the paper is organized as follows. Section ^] contains our results on the ordering in the 2D nn and 
nnn Ising model at fl = and ^ 0. The results on the spin-SAF transition in the coupled model are presented in 
Section rrrn The final Section HVI presents the summary and discussion. 



where the bold variables denote lattice vectors, the first [second] sum includes only nearest neighbors (nn) [next- 
nearest neighbors (nnn)] of the lattice, respectively. Spins along the sides of an elementary plaquette interact via the 



II. 2D NEAREST- AND NEXT-NEAREST-NEIGHBOR ISING MODEL 



We consider the 2D Ising Model on a square lattice with the Hamiltonian 




(1) 
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nn coupling Jn, while spins along plaquette's diagonals interact via the nnn couplings Jki = J1.2 (see Fig. 0. The 
way we defined the Hamiltonian corresponds to antiferromagnetic couplings for Jj > and ferromagnetic for < 0. 

A. Ground state phases 

There is no exact solution of the model Its possible ordered phases and critical properties have been studied 
within various approaches for equal diagonal couplings Ji — J2 (see f?] for a review and references on the original 
literature). We will consider arbitrary Ising couplings (Ji, J2, Jn), so the model can be either frustrated or not 
(see footnote^. The ground state phase diagram can be found from energy arguments, as was first done by Fan and 
Wu for Ji = J2 jSi]- (Their phase diagram is shown in Fig. From direct counting of the ground state energies of 
possible spin arrangements we construct the phase diagram for Ji ^ J2. Along with the phases found by Fan and 
Wu- ferromagnetic (FM), antiferromagnetic (AF), and super-antiferromagnetic (SAF)- there is a fourth phase which 
can occur if Ji and J2 have opposite signs. The name of the SAF phase comes from viewing it as two superimposed 
antiferromagnetic lattices (one lattice of circled sites and another of squared sites in Fig.^. In the SAF state there 
are two frustrated bonds Jn per plaquette and its energy is four-fold degenerate, since each of the superimposed 
lattices can be flipped independently. In addition to the two SAF states with alternating ferromagnetic order along 
the horizontal chains (one of these is shown in Fig.^, there are two states with the vertical ferromagnetic order. 

The new fourth phase shown in Fig. ^ can be viewed as two superimposed lattices each of which is ordered 
ferromagnetically along one side (e.g., J2 < 0) and antiferromagnetically along the other (e.g., Ji > 0). So we will 
call it super-ferro-antiferromagnetic (SFAF) 9]. The SFAF state also has two frustrated plaquette's bonds and four- 
fold degeneracy. The ordering pattern shown in Fig. ^ changes only by a lattice spacing shift over flipping of the 
sublattices. The direction of the ferromagnetic order is determined by the ferromagnetic diagonal. 

The ground state phases in the space {Ji,J2,Jn) are shown in Fig. |21 In order to facilitate perception of this 
picture, we also present in Figs. I3I4I several plane projections of the 3D Fig.[2| In the first quadrant (Ji, J2 > 0) in 
the region Ji + J2 > \Jn\ lying between two frustration planes FP 

Ji + J2 = |Jn|: FP (2) 

the ground state of the model is SAF. A continuous transition from the paramagnetic (PM) to the SAF phase occurs 
at some critical temperature Tc > 0. From the arguments known for the case Ji = J2 UJ (see also ll3| for a more 
general symmetry analysis of the Ginzburg-Landau functional) this transition is non-universal: the critical indices 
continuously depend on the couplings Jn, Ji, J2. 

In the region Ji > Jn of the fourth quadrant ( Ji > 0, J2 < 0) lying between the other pair of frustration planes 
FP' 

Ji = |Jn| : FP' (3) 

the ground state is SFAF. The same arguments 0^,1^ suggest non-universality of the PM SFAF transition. 

In the regions lying above (beneath) the frustration planes FP and FP', and above (beneath) the basal plane Jn = 
in the third quadrant, the ground state is a usual AF (FM), respectively. The transition PM AF (FM) belongs to 
the 2D Ising universality class. The second quadrant Ji < 0, J2 > 0, not shown in Fig. |21 is obtained by a reflection 
over the plane Ji = J2- In the AF (FM) state the number of frustrated (diagonal) bonds per plaquette is two in the 
first quadrant, one in the second and the fourth, and zero in the third. 

Transitions at finite temperature should be absent on the frustration planes FP/FP' where the model is highly 
degenerate. We are not aware of studies of the ground state in these cases and cannot say at the moment whether the 
system possesses some kind of a long-range order at zero temperature or not, except a rather trivial line Jn = Ji = 
0, J2 < of the FP' planes crossing where the model becomes a set of decoupled Ising chains, and four special lines 
on the FP planes where it becomes the exactly-solvable Isotropic Triangular Ising (ITI) model. The latter case will 
be discussed momentarily. 

B. Exactly-solvable limits 

In the 3D space (Ji, J2, Jn) beside the frustration planes FP and FP', there are three special planes where one of 
the couplings is zero. On these planes the model (Q) reduces to the exactly solvable cases. 
Let us start with the upper part ( Jn > 0) of the SAF region 



Jn < Ji -I- J2 



(4) 
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On the SFAF-SAF boundary J2 = the model is equivalent to the anisotropic Ising model on a triangular lattice 
(ATI), for which exact results are available |14L llsL Ha . [l7l|. The antiferromagnetic (Jn, Ji > 0) ATI model with one 
strong bond Ji > Ju is disordered at any non-zero temperature 0|. It orders only at T = 0, i.e., it is Quantum 
Critical (QC). The highly degenerate ground state (however with a vanishing zero-temperature entropy per site) can 
be viewed as a 2D array of antiferromagnetically ordered (along the strong bond Ji) correlated chains. The oscilla ting 
(with a period of four lattice spacings) power-law decay of the spin-spin correlation function along Jn-directions |l7| 
indicate on the preference of the ferromagnetic order along the "missing" diagonal J2. This resembles the SFAF 
state, however any couple of adjacent Ji-chains is uncorrelated. We label this state occurring on two sectors of the 
Ji(or J2) = planes as ATI (QC) on the phase diagram (Fig. [5J. Since the critical behavior of the ATI model with 
two equal weak ferromagnetic bonds | Jn | < Ji is equivalent to the totally antiferromagnetic ATI model , the ATI 
(QC) state smoothly continues into the lower ( Jq < 0) part of the SFAF-SAF boundary. 

The sectors J2 = 0, Ji > (and 1^2) above the FP Ji -I- J2 = Jn correspond to the antiferromagnetic ATI 
model with one weak bond Ji < Jn. It is known [T^ to have only two phases and to order at finite temperature. 
Tc(PM AF) as a function of couplings is also known exactly. The sectors J2 = 0, Ji < 0, Jn > (and 1^2) 
correspond to the antiferromagnetic ATI model which is even not frustrated, and Tc(PM — > AF) > at any Jn > 0. 
The PM AF transition in the ATI model belongs to the 2D Ising class. So, except the SFAF-SAF boundary, the 
ordered phase on the exactly-solvable "triangulation" planes Ji(or J2) = is the same as the AF ground state in the 
interior in this region of the phase diagram. 

The situation on the "triangulation" planes in the lower part ( Jn < 0) of the phase diagram is exactly analogous 
to the upper part, with an obvious replacement AF FM. 

Note that the ground states change on the lines where the triangulation and frustration planes cross. To put it 
differently, these are the lines of quantum phase transitions. The AF (FM) phase disappears in the limit Ji 
\Ju \ — (J2 = 0) . Also, the zero-temperature AF in-chain order [ATI (QC)] described above disappears in the limit 
Ji ^ |Jn| + (J2 = 0) as well. The ITI model Ji = Jn is disordered at any non-zero temperature (indicated as ITI 
(QC) in Fig. EJ. Its ground state, albeit having finite entropy per site, possesses periodical (with a period of three 
lattice spacings) long-range order [l^ . 

The basal plane Jn = in Fig. |21 corresponds to the case when Hamiltonian (QJ represents two decoupled identical 
nn Ising models residing on two superimposed lattices (shown by circles and squares in Fig.^. Diagonal couplings 
J1.2 are the nn couplings of these Ising models. This is the only exactly-solvable limit (labelled by 2 x 2DI in Fig. 
121) within the SFAF (or SAF) region of the phase diagram. In this limit the PM SAF (or SFAF) phase transition 
enters into the 2D Ising universality class. 

C. Incommensurate (floating) phase 

So far we have discussed the ground-state phases of the model and the critical behavior on the boundaries of these 
phases with the disordered phase, as well as on the special planes (lines). However, there is also a possibility that 
ordering into the ground-state phases of Fig. [21 happens not necessarily from the PM phase, but from some other 
one(s) occurring at non-zero temperature. A very simple analysis indicates that this indeed can take place in our 
model. Fourier-transforming the Hamiltonian Q we obtain (we set the lattice spacing to unity) 

= 5]j(q)T-(q)T-(~q) , (5) 
q 

J(q) = JnicOS q^j: + cos Qy) + JiCOS^q^: - Qy) + J2C0s{q3: + Qy) 

where q runs within the first Brillouin zone \qx.y\ < tt. At mean-field level, a minimum of J(q) in q-space defines 
the wave-vector qo of the critical freezing mode r^(qo), i.e. the order parameter (T^) cx cos(qom + ip) below a 
certain critical temperature Tc- In different regions (Ji, J2, Jn) of the ground-state phase diagram (Fig. (21) we find 
minimum at q''''* — (0,0)/(7r,7r) giving the FM/AF order parameter and two minima qf'^2'' = (""j 0)/(0, tt) giving two 
components of the SAF order parameter. (The latter represent two possible ordering patterns of the SAF phase and 
via some transformation can be related to magnetizations of the superimposed sublattices ) . The important point 
is that the positions of the commensurate (C) extrema q" (jj =F,A,SAF) of J(q) do not depend on couplings. 
There are however two other pairs of extrema ztq"'" which exist if 

|Jn|<2|Ji| and/or 2|J2| (6) 

These extrema lie on the diagonal of the Brillouin zone. q*'° are generically incommensurate (IC) and depend on 
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couplings as 

Qx ^ <ly^ arccos ( - ^) (7) 

= = arccos ( - ^) (8) 

We show the positions of aU extrema of J(q) within the BriUouin zone in Fig. |S1 In the SFAF ground state region 
(e.g., in the fourth quadrant Ji > 0, J2 < shown in Fig. |2J) the pair of extrema ±q" (|HJl gives global minima of J(q) 
(the other pair of solutions © ±q* when exists, corresponds to its maxima), and \q^\ — \qy \ > 7r/2. As we see, two IC 
modes ±q" could give components of the SFAF ground state order parameter's wave vectors iq^''*'' = ±{tt/2, —tt/2) 
only if Jo — 0. (In the second quadrant of the SFAF region when Ji < 0, J2 > 0, the vectors q'' and q° exchange their 
roles. Because of the Ji J2 symmetry, in the following we will always discuss the fourth quadrant for concreteness.) 

The locus of the IC global minima does not coincide with the SFAF ground state region, but overlaps with the 
neighboring FM, AF, and SAF phases. The minimum J(±q°) is located between two planes | Jn| = 2Ji in the fourth 
quadrant, and in two regions of the half of the first quadrant (J2 < Ji): (i) between Jn = +2Ji and Jn = +2\/ J1J2] 
(ii) + i-^ — . The regions of the IC minima in the other half (J2 > Ji) of the first quadrant (as in the second quadrant) 
are obtained from the described above by Ji ^ J2 , q° 1-^ q" . On the plane phase diagram shown in Fig. 0] this locus 
is restricted by the lines y = 1/2, x = 1/2, y = l/ix, shown by the black dashed lines. 

So we can conclude that at finite temperature the model possesses an IC phase and there is an IC-C phase transition 
where the IC wave vector q" locks into one of the (commensurate) ground-state phase vectors. As in 2D the IC phase 
has only an algebraic long-range order, it is called floating The origin of the IC floating phase in our model is 
frustration (competing interactions). Such phase is well known from another example of frustrated Ising model, i.e., 
the ANNNI model which was intensively studied in the past P, ITol ITU . In that model the floating phase locks into 
the antiphase which has the wave vector q — (0,7r/2). (The antiphase is analogue of our SFAF phase.) The ANNNI 
model also provides an example showing that the mean-field (minimization) analysis does not work well in defining 
boundaries between the floating and commensurate phases in 2D, and the extend of the IC phase is less than the 
mean field suggests ^. 

We will not attempt to locate exactly the phase boundaries at finite temperature in this study. Following Domany 
et al |l3j in classification of the ordered phases by commensurability p, i.e., the ratio of superstructure's period and 
lattice spacing along a given direction, we can label the phases as follows: F > 1x1; AF i-^- 2 x 2; SAF 1-^ 1x2 
(or 2x1); SFAF 4x4"^. From mapping of the 2D IC-C phase transition to the Kosterlitz-Thouless problem, 
it is established that there is no such (continuous) transition for commensurate phases with small < 8 0, . 
From this result with the proviso of continuity of phase transition(s) in the model, we conclude that the IC floating 
phase cannot "spill" beyond the SFAF (j) — 4) ground-state region of the phase diagram, even if a naive (mean-fleld) 
analysis suggests that within the F, AF, and SAF regions there are some parts where it could be possible (see Fig. 
2J. The only high-temperature phase the latter three regions have a common border, is the disordered FM. 

Combining this with the known exact critical properties of the model on the special planes discussed above, we end 
up with the qualitative flnite-temperature phase diagram shown in Fig. El Since on the plane Jn = the model ^ is 
just two decoupled Ising lattices, the floating phase must be absent. Mean-fleld arguments suggest that the floating 
phase disappears exactly at Jn = giving rise to a Lifshitz point (L). Note that the floating phase does not appear if 
the diagonal couplings are equal even at the mean-field level (see Fig.^J, which agrees with known more sophisticated 
analyses of this case |l|. 

D. The model in transverse field 

Now we turn to the analysis of the nn and nnn Ising Hamiltonian i? in the presence of a transverse field. The 
total Hamiltonian of the Ising Model in Transverse Field (IMTF) reads 

HiMTF = H - (9) 
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The well-studied ANNNI phase diagram with the FM and antiphase ground states Q, llll is an analogue of the lower part of the 
fourth quadrant (FM-SFAF) of our diagram in Fig. ^ 

According to this notation, the antiphase of the 2D ANNNI model is (1 X 4). For the same reasons we give for our case, the floating 
phase of that model exists only within the antiphase ground state region. 
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The Ising operators are normalized to satisfy the spin algebra 

[TC,T[]^^6,,eo.p,T^ (10) 

There is no exact solution of the transverse 2D Ising model even for the case of nn couplings only ( Ji = J2 = 0). The 
ground state phase diagram of the Hamiltonian iPJ can be analyzed from the known mapping of the d-dimensional 
IMTF at zero temperature onto the (d + l)-dimensional Ising model at a given (non-zero) "temperature" 3]. In our 
case the 2D nn and nnn IMTF maps onto the 3D Ising model comprised of the 2D layers coupled ferromagnetically 
in the third (Trotter) direction with the coupling Jt (x — IncothO < 0. For such (2 + l)-dimcnsional model a mean- 
field analysis gives a qualitatively correct diagram of the ground state phases of the 2D IMTF ^ . The new coupling Jx 
does not bring any additional frustration to the 2D nn and nnn model. Analysis of Jsi^Qx, Cly, Qt) = Jt cosqT + J{qx, %) 
where J{qx, %) is given by (jsj, shows that Jt does not modify the domains of the global minima in the ( Ji, J2, Jn)- 
space, adding only a trivial qr = third component to the two-dimensional vectors discussed above (cf. Fig. 
0) The temperature phase diagram of the 3D Ising model with the spectrum J3{qx,qy,qT) (if we label the phases 
according to the in-plane ordering pattern defined by the 2D vectors q^) looks similar to the one shown in Fig. El with 
one very important distinction: the abovementioned argument related to phase's commensurability p does not apply 
in 3D, so the IC region is not restricted to lie above the SFAF phase, but can spill into the neighboring regions of the 
{•h, J2, Jn)-space. From the mean-field arguments, the IC region is given by the locus of the IC minima J(±q°/'') 
defined in the previous section. So, the IC phase instead of being locked between the special planes FP' and J2 — 
as shown in Fig. [S] (a) and (b), respectively, can spread up to the locus boundaries shown by the crosses. From the 
equivalence between the zero-temperature d-dimensional IMTF (quantum Ising) and the {d+ l)-dimensional classical 
Ising model, we infer that the ground-state phase diagram of the former on the plane {^/ J, Jj) should have the same 
structure as the described above (T, Jj) diagram of the latter. So we expect the transverse field to generate the IC 
ground-state phase not only in the SFAF region of the the Ising coupling space, but also in the neighboring parts of 
the F, AF, and SAF regions. From minimization arguments the latter are restricted by the dashed lines on the plane 
diagram in Fig. ^ and by the crosses in Fig.El 

From analogies with the ANNNI model, we rather expect this "IC region" to be filled with infinitely many com- 
mensurate phases with different p but detailed analysis of this question, as well as the full finite-temperature 
phase diagram of the transverse model @ need a separate study. 

In the rest of the paper we will be particularly interested in the SAF region of the coupling space, and restrict 
ourselves to the couplings 

j2<4JiJ2, (Ji,J2)>0 (11) 

According to Fig. ^ it means that we choose the couplings to lie above the hyperbole y = 1/Ax. The first condition 
in (|ll|l ensures that the couplings lie in the region where qf*2'' = (0, 7r)/(7r, 0) provide a global minimum of model's 
spectrum, so the phase with the IC solutions q"/* does not intervene. The second in the above conditions stipulates 
that even if Jo = we stay away from the planes where our model becomes the triangular Ising. A transverse field 
can can generate exotic temperature phases in that model. Such phases were found 0, Q in the particular case of the 
isotropic (antiferromagnetic) transverse triangular Ising model"^. 

From mapping between the quantum and classical Ising models we conclude that at zero temperature our IMTF 
with the couplings satisfying Hll() possesses a single Quantum Critical Point (QCP) which separates the SAF and 
PM ground state phases. The mean field predicts a two-phase PM/SAF diagram. The critical temperature Tc of the 
second-order PM-SAF phase transition evolves smoothly from the QCP Tc{J/^a) = to the asymptotic limit Tc(oo) 
of the classical model (see dashed curve in Fig.|Sl where g = (Ji -I- J2)/fl). It is also known that the mean field gives 
a qualitatively correct phase diagram for the IMTF when d > 2 i3|. Thus we argue that the mean-field result shown 
by the dashed curve in Fig.|Slrfoes represent the phase diagram of the IMTF H1I9I11|I . while its quantitative aspects, 
e.g., the exact value of the QCP, should be corrected via more accurate treatments. 



By analogy with the triangular Ising case of Refs.^l^, we expect the transverse field to bring about new exotic phases near the 
frustration planes FP, FP', where the ground state of our model is also infinitely degenerate. 
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III. COUPLED SPIN-PSEUDOSPIN MODEL 



A. Hamiltonian 



Now we turn to the analysis of the IMTF coupled to the quantum spins (S) residing on the same sites of 

the lattice as the Ising spins do. The latter we will call pseudospins from now on. Such coupled spin-pseudospin 
(or spin-orbital) models emerge in various contexts, most notably the Jahn- Teller transition-metal compounds [l9l| or 
many kinds of low-dimensional quantum magnets. For a recent short overview and more references, see 0. 

Models of the type we study in the present work appear in analyses of the quarter-filled ladder compound NaV205. 
The Hamiltonian of this material can be mapped onto a spin-pseudospin model with spins and pseudospins residing 
on the same rung of a ladder |23,|^|23. The ladders form a 2D lattice. The long-range pseudospin order (T^) ^ 
represents physically the charge disproportionation between left/right sites on a rung below 

This system was analyzed on the effective triangular lattice shown in Fig. [7^ by solid lines |20ll2^ . However, in the 
case of NaV2 05 the Ising couplings generated by Coulomb repulsion are antiferromagnetic and Ji > Jq. Since the 
triangular Ising model with one strong side is disordered 0, one needs an extra coupling ( J2-diagonal) to stabilize 
the observed SAF long-range order. In our earlier study we explicitly took into account Ji, while the other diagonal 
J2 was effectively generated via the spin-pseudospin coupling^. 

In this work we take into account the Ising couplings Jn, Ji, J2 between neighboring sites of the effective lattice, 
as shown in Fig. |7^. Then such effective lattice can be mapped onto the square lattice shown in Fig. [TJa with the nn 
and nnn Ising couplings. For NaV205 all Jj > 0, so the model is frustrated. As follows from geometry of the original 
NaV2 05 lattice, J2 is the weak diagonal and Ji is the largest coupling: 



J2 < Jn < Ji 



(12) 



We assume Jn to be small enough not only to lie beneath the frustration plane but to satisfy a more stringent 
condition Then according to the above analysis, the IMTF with these couplings has a two-phase (PE-SAF) 

diagram with a QCP^. In the literature on NaV205 its charge order is called the "zig-zag phase", what characterizes 
the antiferroelectric order in a single ladder only. In fact, the two-dimensional long-range charge order in NaV205 is 
SAF. For a detailed explanation of this point, including interpretation of the experimental crystallographic data on 
the charge order |2^ in terms of the Ising pseudospins, see 24]. 

In the following we will work with dimensionless quantities: Hamiltonians Ti. = H/Q, temperature T —* T/Q and 
Ising couplings g^ = Jj/ri. With site labelling shown in Fig. [7| the IMTF Hamiltonian is: 



1 



9n 



\ mn 



l,n+l 



rrx rrx 
mn m+l,7i 



91%. 



rj-X 

mn m.n-\-l 



For the decoupled spin sector of the total Hamiltonian we take into account only the strongest coupling between 
spins on the nn rungs of a ladder. In terms of the effective lattice (cf. Fig. ^ this translates into a set of decoupled 
Heisenberg chains with the usual antiferromagnetic spin exchange J. These parallel chains are oriented along the 
Ji-diagonals. 

As we infer from our previous work on a simpler version of the IMTF Hamiltonian , there are two spin-pseudospin 
interaction terms resulting in two qualitatively distinct aspects of model's criticality: the inter-ladder spin-pseudospin 
interaction oc e, and the in-ladder spin-pseudospin interaction oc A. The former, in terms of the equivalent square 
lattice, linear over difference of the charge displacement operators on the nnn sites along the weak J2-diagonal, 
is responsible for the simultaneous appearance of the SAF order and the spin gap, as well as for the destruction of 
the IMTF QCP^. The latter, quadratic over the nnn charge operators along the Ji-diagonal, is responsible for the 
re-entrance. With all these terms the total effective Hamiltonian reads: 



"H — TiiMTF 



'm,n-\-l 



[J- 



^ ran m,n+l 



(14) 



^ A linear coupling oc <f(T\ T ^) with some Gaussian mode i/s results in an effective anti-/ferrro-magnetic interaction between Ti and 
Ti. Such terms can create or renormalize the couplings of the Ising effective Hamiltonian. In his context the Ising model with two nnn 
couplings of different signs is less academic than it might appear. 

® Since in applications to NaV2 05 the Ising pseudospins T represent charge displacements, the appropriate names for the phases are 
"paraelectric" (PE) and "super-antiferrolelectric" . We keep the same abbreviation SAF for the latter. 

This type of coupling is allo wed by the symmetry of the original NaV205 lattice 0|. Numerical estimates of e from a microscopic 
Hamiltonian are given in |25(. Note also that e effectively couples the spin chains. 
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The dimensionless couplings J, A, e in the Hamiltonian H14() are positive, and the spin operators satisfy the same 
algebra H1U|) as the pseudospins (while S and T commute). The sums above run through 1 < m < A4 and 1 < n < Af. 
For brevity we will use the notation Dmn = Smn^'m,n+i- In this study we consider the model with the XY spin 
sector: 

Dmn = *5'„„S'„ „^]^ + 'S'^jjS'^j (15) 

The range of couplings under consideration will be restricted to 

iJA)<gi, e<max{J,A} (16) 



B. Spin-SAF phase transition 

We treat the Hamiltonian p4|l following conventional wisdom of molecular- field approximations (MFA) . In the 
present version of MFA the pseudospins are decoupled and averaged with the density matrix cx exp{— phmn'^ mn) , 
where hm„ is the Weiss (molecular) field, while the spin sector is treated exactly via a Jordan- Wigner transformation. 
The details are presented in 6]. Similar to the pure IMTF with couplings Hll|l we assume the possibility of the SAF 
order in the coupled model (|14|l . So we take the following Ansatze for the Ising pseudospin averages (i.e., the charge 
ordering parameters in terms of the real physical quantities) 

(^™„> = ™. (17) 
- (-l)™+"m. (18) 

It is easy to see from the Hamiltonian ()14|l that ansatz H18|l creates a dimerization in the spin sector, therefore a 
natural assumption for the dimerization operator average is 

{D„,,,)=-[t+i-ir+"S] (19) 

With the new coupling 

9 = 91 + 92 (20) 

the molecular-field equations and results derived in for the case 92 = 9n = (i-e., 9 = 91) can be applied here. 
Some of them we reproduce in this paper in order to make it more self-contained, and for the use in what follows as 
well. 

The average quantities are determined by the system of four coupled equations 

ll + 2Xtrn, 131) 

niz = tanh — (21a) 

m^g + 2er] B\) 

rrix = — tanh— (21b) 

t = - dip— ^ tanh p^{ip) = -tniA,P) (21c) 

A sin^g? _ A - 
77 = / dip— ^ tanh pa^) = Vn{A,(3) (21d) 



where f) — yTif+Tif is the absolute value of the Ising molecular field 

h,_ = 1 + 2Xtmz (22a) 

hx = 9mx + 2eS (22b) 

The other auxiliary parameters are defined as follows: 

S = mxTj (23a) 

^(^) EE \J cos2 93 -f A2 sin^ ip (23b) 

P = ^{J + XmD (23d) 
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At some critical temperature Tc the coupled model undergoes the phase transition. It is of the second kind, with the 
thermodynamic behavior of the physical quantities as the Landau theory of phase transitions predicts |^ . With the 
spin-pseudospin (-charge) coupling e present, the SAF charge order nix 7^ and the spin gap Asg = 2emx appear 
simultaneously below Tc. By analogy with the spin-Peierls transition, when the Peierls phonon instability (freezing) 
creates the spin gap, it is natural to call this type of transition the spin- super- antiferroelectric (spin-SAF) transition. 

It is worth to point out an important property of the Hamiltonian (|14|l : in the other domains of Ising couplings 
(not considered in our analysis of the coupled model) where the Ising sector of (|14f) can order into, e.g., FM, AF, or 
SFAF phase, the dimerization (gap) in the spin sector does not occur. 

The behavior of Tc{g) in the coupled model H14|l shows two new striking features comparatively to the pure IMTF: 
re-entrance and destruction of the QCP '6'|. In the absence of the spin-charge coupling e, the model (|14l) has a QCP 
at 

3,^2(1 + ^), (24) 

where Tc vanishes (see Fig.lSJ. A renormalizes the QCP comparatively to the pure IMTF value g ~ 2. The coupling e, 
responsible for the spin gap generation also destroys the QCP, resulting in the exponential behavior of Tc in the region 
g gx, where the model would have been disordered at any temperature if e = 0. This constitutes an important 
feedback from the spins on the charge degrees of freedom, allowing the very possibility of the model to order at all. 
Approximate analytical solutions for Tc{g) in the regimes of strong Ising couplings and the BCS- are: |^ 



« <^ - - (25) 

4J5 



Mexp[ - ^(ffA - 5)] , BCS regime 



where A = — ~ 1.6685, 7 ~ 0.5772 is Euler's constant, and 

J^J+^ (26) 
4 

The boundary where the low-temperature BCS regime sets in and the related formulas are applicable, is given 
approximately by the condition 

4e2 

BCS regime : g < gx ^ ^ (27) 

ttJ 

The BCS regime has many analogies with the standard theory of superconductivity, apart from the exponential 
dependence of Tc on couplings. In particular, several physical quantities (order parameter, BCS ratio, specific heat 
jump) manifest certain "universal" behavior near Tc, similar to that known from the BCS theory |6|. 

Another particularity of Tc{g) found earlier |^ from the numerical solution of Eas. l|21|) . is re-entrance in the 
intermediate regime g ^ gx- The re-entrance occurs in the coupled model with the QCP (e = 0), while when e 7^ 
the critical temperature can even manifest a double re-entrant behavior before it reaches the BCS regime (see Fig.|SJ|. 
A detailed analysis of the coupled model in the regime of re-entrance was not done previously. We will address this 
in the next subsection, mainly analytically, in order to get more insight on the underlying physics and, in particular, 
to establish conditions when the re-entrance can occur. 



C. Re-entrance 



Let us first reproduce some earlier formulas [g for reader's convenience. At T < Tc one equation from the pair 
()21b.b') can be written in a form 



At T = Tc we have Eqs.lEHi) as 



1 2A , , 

m-'^g+ ,, r/„ i„, T<Tc (28) 



1 , 0c , 2Amr 
= -tanh^(l + U„) , (29) 

Z Z TT 
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and parameters t„,r7„ are given by Eas. 1)211 . d) with A = 0. The latter two functions have the following expansions: 
i 

tn{Q,x)^{ (30) 



and 



r ^x{l - ^x^) + 0{x^), X < 1 
Vn(.0,x)^\ (31) 
[ lnAa; + fg^+0(^), x>l 



1. Case e = 0; re-entrance with QCP 

As one can easily see from Eas. (|28l29|l there is no re-entrance when A = 0. This is a well-known fact for the 
pure IMTF, as on the mean-field level, as well as beyond MFA 0, |2^. To study the re-entrance analytically and in 
particular, to establish whether there is some minimal value of A when it appears, we should distinguish between two 
asymptotic regimes of the mean- field equations. Let us first consider the regime (it can occur only if J < 1) when 

\j <Tc<\ (32) 

(In all regimes of couplings the re-entrance occurs at Tc < \)- By carrying out the leading-term expansions of the 
functions in Eas. (|28l29l) we obtain 

g = 2 + 4e-V^^ + ^ (33) 

for a single-valued function g{Tc). The non-monotonic (i.e., re-entrant) behavior of Tc{g) is related to the existence 
of an extremum of g{Tc). The coupling gmin which defines the minimal value of g for the order in being possible 
(in the pure IMTF with A = this was the QCP), and in the same time the left border of the re-entrant region 

.gmin < g < g\, (34) 
is defined from the minimmn of the function g(Tc) (|33|l . This point corresponds to the critical temperature 

T, — In^^ Ko, Ko = — ^ (35) 
AJ 



for which 

4 



> + lnKo) (36) 



The consistency of the solution H35|) with l|32|l implies the condition 

2<lnKo<^, (37) 
J 

The other regime corresponds to the case when 

Tc < min{l/2, J/2} (38) 
Proceeding in the same way as above, we obtain for g{Tc) in this case: 

g^g,+ 4e-^^/2^^ - (39) 

Let us point out that the conditions (|32I38|I determine two different regimes (x < 1 or x > 1) of the asymptotics 
(|30I31|I we apply in order to obtain g{Tc) as (|33|) or (|39|l . So if ,/ < 1 there are regions of Tc where condition H32|) 
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is satisfied, then the approximation (|33|) apphes. However at sufficiently low temperatures {Tc < J/2) we inevitably 
enter the other regime (I38|) where the asymptotics (|30l31f) change (a; < 1 i-^ x > 1), and the function g{Tc) crosses 
over from H33|l to H39|l . If, on the contrary, J is large, then condition H32|) never applies, and the approximation H39|) 
describes the whole region Tc < 1/2. 

Extrema T, of the function g{Tc) (|39|l are determined by the transcendental equation 

g-<,./2T, ^ (40) 

This equation always has a trivial solution = corresponding to the (local) maximum of g{Tc). This is the QCP, 
and the curve Tc{g) approaches the QCP normally to the abscissa (see Fig.|SJ). Two non-trivial solutions of 14UI) exist 
if the couplings satisfy the condition 

J>CiX-^gx, (41) 

where 



^(0.0.3121. (42) 



There is only one solution within the validity region of the approximation H39() . and it corresponds to the minimum 
of g{Tc). If the couplings satisfy (|41|l then 

u = ln3+ ilnoo > 1 (43) 
24J2 



and the minimum can be found analytically as 



T.^'-^C--'^). (44) 



For the left border of the re-entrant region we obtain 

.,™^..-^T.^ + ^rJ (45) 

The above equations agree well with the numerical solutions of the MFA H21|l at different values of couplings (from 
comparison of the asymptotics (|33l I39|l and the numerical curves at various couplings and temperatures we found 
the deviations ^ 5% at most). More importantly, the analytical results of this subsection allows us to understand in 
details the interplay of the scales provided by model's couplings and the temperature, resulting in the re-entrance. 
Let us explain this on the example of two characteristic numerical curves shown in Fig. |H| 

The curve shown for J = 0, A = 1, e = ( J = 0.25) corresponds to the case of small J. At Tc > J/2 ~ 0.125 it is 
well described by the equations for the regime H32II . Its re-entrant behavior and, in particular, the minimum gmin is 
due to the last term on the r.h.s. of (|33l) . At lower temperatures Tc < 0.125 the asymptotics (|33f) is not applicable, 
the curve is described by H39|l . Note that since CiX^gx « 0.8229, the condition (|41l) for the minimum is broken, and 
the asymptotics (|39|l describes the featureless low-temperature evolution of this curve towards the maximum at the 

The second curve in Fig. |Slwith J — 0.75, A = l,e = (J = 1) corresponds to the case of large J. The whole 
re-entrant region (Tc < 0.4), including the position of the minimum (jmin {CiX^gx ~ 0.8229) is described by the 
interplay of the last two terms on the r.h.s. of Ea. (|39|l . Comparison of Eas. H3t)l45|l allows also to understand a more 
pronounced re-entrant behavior for the case of smaller J. 

As we see from our analysis of Eas. (|33l39|l in the both regimes (|32I38|I . the the re-entrant behavior on the phase 
diagram occurs at any X 0. 

2. Case e 7^ 0; double re-entrance, no QCP 
The absence of re-entrance at A = can be proven rigorously. Indeed, combining Eas. H28l29|l we obtain the equation 

4e2 J, l + 2m. 



m. 



-=.+ ^.„(0,2ln^) (46) 
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which has one and only one solution niz G [0, 1/2] for a given value of g. This solution in its turn provides a unique 
value of Tc via Ea. H29|) . thus no re-entrance. 

At A 7^ continuous evolution of Tc{g) between the regimes of strong Ising coupling and BCS- [cf. Ea. (|25|l ] can 
occur either with a double re-entrance [i.e., with one minimum and one maximum of g{Tc)] within the re-entrant 
region 

(47) 

or without re-entrance. In the latter case the function Tc{g) [or g{Tc)] has only an inflexion point (see Fig.|Hl). 
Following the analysis given in the previous subsection, we obtain for the case of small J in the regime ()32fl 

5 = 2 + 4e-V^= + ^:^^ (48) 

Again, the re-entrant behavior is conditioned by the existence of a minimum of g{Tc). It exists if, at least 

e < ^ VaJ w 0.7071 VxJ (49) 
V2 

The unique minimum of g{Tc) l|48(l defines the left border of the re-entrance region gmin and corresponds to the critical 
temperature T, given by Eas. H35l36|l with and 



A J - 2e2 



The consistency imposes the constraint analogous to (I37II . more stringent than the "minimal requirement" H49|l . 
The other regime ()38f) is described by the approximation 

g = + 4,-../2T. _ ^ J + (51) 

As we have explained in the previous subsection for the case e = 0, the asymptotics H48|l is applicable only for small 
J at the intermediate temperatures (|32|l . while (|51|l can be applied at arbitrary low temperatures, including the BCS 
region. The latter, given by the exponential dependence in H25() . can be recovered if we retain only the first and last 
(leading) terms on the r.h.s. of Ea. (|51|l . 
In the regime ()38|l the re-entrance occurs if the equation for extrema of H51|l 



{XJ + e^)T^ - (52) 



has non-trivial solutions. Note in making comparison of Ea. H52|) to its counterpart (|40|l at e = 0, that coupling e 
destroys the QCP 0, as immediately seen from (|51|l . So the trivial solution = of Ea. H52l) corresponds to the 
unphysical singularity of g. As follows from H52I38|) . non-trivial solutions are possible if, at least 

£ < V , =v6J?^ 0.8357\/aJ (53) 

If this condition is satisfied, the transcendental equation (|52|l has at least one solution, corresponding to maximum of 
g{Tc). To leading order in e, it occurs at the temperature 



and the coupling 



4^2 AttVAJ 

g^..^g.~^\n^^ (55) 



If the couplings meet both the conditions (|41|l and 



e<C2ffA(j)', (56) 
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where 

then a second solution of H52|) (T"), corresponding to minimum of g{T^) exists. This minimum, located between 

t: < t: < ^ (58) 

is given approximately by Eg . (|44|l where u in H43() is modified by Aq i— > a, and 

24 7-^ 

7 (59) 



For validity of expansion 144|) we assume u > 1. 

The analytical results of this subsection allows us to describe the behavior of the coupled model at e 7^ 0, following 
from the MFA equations l|2U, both qualitatively and quantitatively. In Fig. |H1 two counterparts (e = 0.1) of the 
numerical curves discussed in the previous subsection are shown. For this case of small e, re-entrance is possible, 
according to conditions (|49I53|1 . The re-entrant behavior at the temperatures Tc > Tl is not modified essentially by 
the presence of new coupling e comparatively to the case e = 0, and is in fact controlled by couplings J, A, g. Since 
we have already discussed it in detail for the case £ = 0, we will not dwell on it any more. In the low-temperature 
regime (|51|) coupling e changes drastically the behavior of g{Tc) at Tc ^Tl, creating a maximum at g{Tl) described 
well by the approximation (|54|l and turning g{Tc) away from the QCP towards the BCS region. For the BCS regime 
Ea. (|25|l provides a virtually exact solution. 

As follows from the inequalities H49I53I56|I an increase of e can suppress the re-entrance even at A 7^ 0. The necessary 
conditions H49I53|) for extrema of the two asymptotics (|48I51|) are close, albeit with a rather small mismatch of the 
coefficient. Conditions for re-entrance are more stringent, since they require the consistency between the solutions for 

extrema and the validity ranges of the appropriate asymptotics. The (overrated) critical value £0 ~ 0.7\/xj gives a 
good simple estimate for the boundary where the re-entrance disappears from the whole curve g{Tc), whether J > 1 
or J < 1. An example of the curve Tc{g) without re-entrance is shown in Fig. |S1 

To summarize our analysis of the re-entrance for the cases £ and £ 7^ 0: it reveals the robustness of this 
phenomenon in the coupled model H14I) and its underlying mechanism, namely, competition between different scales 
defined by the couplings J, A, e, g and the temperature. These competing scales (interactions) are not related to the 
Ising frustration which is present in the model as well [(Jn, Ji, J2) > 0], since the latter is not accounted for explicitly 
by our mean- field equations. This competing mechanism for the re-entrance appears to be robust and not being an 
artifact of the MFA. Re-entrant phases due to competing interactions are known also, e.g ., from exact solution of the 
Ising model on the union-jack lattice or from analyses of decorated Ising models pq . 

It is not clear for us at the moment how the proposed re-entrance can be observed. NaV205 is, up-to-date, the 
only known compound with the spin-SAF transition and does not show re-entrance. This is in agreement with our 
estimates for the parameters for the effective Hamiltonian for this compound. They give its g located on the disordered 
side of the (destroyed) QCP, and the re-entrance on the whole curve Tc{g) would be only very weak (i.e., localized 
near g g\), if any. It appears experimentally that, e.g., external pressure cannot modify the in-plane parameters of 
NaV2 05 strongly enough, such that re-entrance would be generated. The variations of the interlayer couplings under 
pressure, on the other hand, generate various types of order (including devil's staircase) with regard to the plane 
stacking, while the in-plane SAF order remains unaffected ^4 



IV. SUMMARY AND DISCUSSION 



We study the 2D Ising model on a square lattice with nearest-neighbor (Jo) and non-equal next-nearest neighbor 
{Ji,2) interactions. The cases of classical and quantum models are considered. 

We find the ground state phase diagram of the classical Ising model at arbitrary Jo , Ji , J2 . Along with the three 
ordered phases -ferromagnetic (FM), antifcrromagnetic (AF), and SAF- known for Ji = J2 0, in a more general 
case Ji 7^ J2 there is a region of the coupling space with the super-ferro-antiferromagnetic (SFAF) [or (4 x 4)] ground 
state phase and an incommensurate (IC) phase at finite temperature, not reported before. The three phases - SAF, 
SFAF, IC - can occur only in the presence of competing interactions (frustrations) on the Ising model's plaquette. 

A particularly interesting conclusion from the analysis of the quantum model's phase boundaries is that transverse 
field {fl) can stabilize the IC ground-state phase (located for = in the region with sign(Ji/J2) = —1) in some 
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parts of the AF and SAF regions of the coupling space where (Ji, J2) > 0, but Ji ^ J2. These regions, along with 
vicinities of the special planes of degeneracy (triangulation and frustration) in coupling space, are good candidates for 
the quantum model to demonstrate a very non-trivial critical behavior. Leaving this for a future work, we hope that 
our findings will inspire additional interest in this model. Taking into account only one simple example of a mapping 
shown in Fig. |7| it is clear the model with Ji ^ J2 is not so exotic. 

We analyze the IMTF coupled to the XY spin chains in the restricted (SAF) region of (Jn, Ji, J2) where the IMTF 
has a simple two-phase (disordered-SAF) diagram with a QCP, similar to that of the transverse nn model. Our interest 
in this model is motivated by the problem of the phase transition in the quarter- filled ladder compound NaV2 05. 
The predictions of the mean-field equations for the critical properties of the coupled spin-pseudospin model do not 
differ essentially from our earlier results for a simpler Hamiltonian 1^. Due to the spin-pseudospin coupling e, the 
QCP of the nn and nnn IMTF is destroyed, and in the whole SAF region (|ll|l of Ising couplings the spin-pseudospin 
model undergoes the spin-SAF transition. We should point out, that albeit the exponential BCS regime H25|) on the 
disordered side of the IMTF QCP g < g\ formally extends up to 5 = 0, decreasing i.e. Ji -I- J2, will eventually 
remove us from the coupling region Hll|l where the SAF pseudospin solution of the MFA equations is applicable. 

We perform a detailed analytical study of re-entrance in the coupled model. In particular, we establish the conditions 
when it can occur. The analytical results not only agree well with the direct numerical calculations in various regimes, 
but allows us to understand the physical mechanism of re-entrance due to interplay of competing interactions in the 
coupled model. 

In this work we gain more insights on the transition in the spin-pseudospin model, and we can sharpen our previous 
statements concerning the applications to NaV205 0. The present analysis of Ising sector allows us to identify the 
2D long-range charge order in that compound as the SAF phase. As follows from known results on the ordering of the 
frustrated 2D Ising model into the SAF phase U, the spin-SAF transition has non-universal coupling-dependent critical 
indices. Experiments indicate rather wide regions of the two-dimensional structural (charge-ordering) fluctuations 
characterized by the critical index (3 « 0.17— 0.19 [s^, 1^112, close to (i = 1/8 of the 2D Ising model. Due to known 
difficulties in extracting critical indices from experimental data, it seems problematic to diagnose the deviations from 
universality caused by < Ju/Ji < 1- (Note that in the limit Jn ^ the PM-SAF transition enters into the 2D 
Ising universality class, and for NaV205, due to its geometry, the ratio Jn/Ji should be small.) 

On the theory side, the critical indices of the PM-SAF transition as functions of couplings in the nn and nnn Ising 
model have been calculated by various methods only at Ji = J2 [l|] . The critical exponents are unknown for the case 
Ji ^ J2, and it appears to be an interesting problem to study. 

Another very interesting issue we addressed recently in a separate study j23|, is the 3D nature of the transition in 
NaV2 05. According to the correlation lengths measurements [33, upon approaching Tc = 34K the 2D 3D crossover 
of the pretransitional structural fluctuations occurs somewhere at T ~ 50K. The model considered in the present work 
deals with a single plane, leaving aside the question of charge ordering along the third (stacking) direction. The phase 
transition in NaV2 05 quadruples the unit cell in the stacking direction, and the recent X-ray experiments, carried 
out deep in the ordered phase [s^l revealed peculiar stacking ordering patterns of the super-antiferroelectrically 
charge-ordered planes. In addition, the pressure can change these patterns and even generate a multitude of higher- 
order commensurate superstructures in the stacking direction (devil's staircase) p9l |. To explain these phenomena 
we proposed a 3D extension of the Ising sector with additional competing couplings between the nearest and next- 
nearest planes j24?|. In the limit Jg — > the Ising sector reduces to two identical interpenetrating decoupled 3D 
ANNNI models. Although inclusion of the competing interlayer couplings accommodates the explanation for the 
observed stacking charge order in the framework of the spin-SAF (in-plane) mechanism of the transition in NaV205, 
a deeper understanding of the critical properties of a very complicated model with the 3D Ising sector warrants a 
further work. 
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FIG. 2: Phase diagram of the ground states (GS) of the model Q. For visuahzation purposes it is drawn within ( — 1,1) x 
(-1, 1) X (-2, 2) parallelepiped. The SAF GS (red) lies between the frustration planes Ji + J2 = | Jn| (bold black). The SFAF 
GS (magenta) lies between the frustration planes Ji = |Jn| (bold black). The AF GS, green (FM GS, blue) lies above (beneath) 
the frustration planes and above (beneath) the basal plane in the third quadrant, respectively. The second quadrant (not shown) 
is obtained by a reflection over Ji — J2 plane. Different sectors shown by hatched and twiggy lines on the exactly-solvable 
planes Jf — are explained in the text. 
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FIG. 3: Plane projections of the 3D diagram shown in Fig. |5| (a): Upper part Jn > 0, view from the top. (b): Lower part 
Jn < 0, view from the top. (c): Compactification of 3D Fig. |5|in the special case J\ = J2 = J' 
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FIG. 4: Plane phase diagram for the ratios of the couplings, corresponding to the 3D Fig. |5]at Jn > 0. The phase boundaries 
(thick solid lines) corresspond to the frustrations planes FP and FP'. The thick dashed lines (black, online) indicate the 
boundaries of the incommensurate global minima locus {y > 1/2, x > 1/2, y < l/ix) discussed in the text. The case Jn < 
can be obtained by the substitutions Ji,2/Jn i—^ Ji,2/\Jn\, AF i— > FM in this figure. 
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FIG. 5: Positions of extrema of J(q) © in the Brillouin zone. Open symbols connected to their bold counterparts by a 
reciprocal lattice vector. Incommensurate extrema iq"'* I|7l8ll (shown for the case Ji > 0, J2 < 0) exist if conditions Q are 
satisfied. 




FIG. 6: Qualitative temperature phase diagram. The IC floating (Fl) phase lies beneath the dashed lines. Crosses indicate the 
borders of the IC minima locus, (a): For fixed Ji > 0, J2 < 0. On the plane Jn ~ the floating phase must be absent. We 
assume that it smoothly disappears at Jn = resulting in a Lifshitz point (L). (b): The same for fixed Ji > 0, jJnl < Ji. 
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(a) (b) 

FIG. 7: (a): 2D effective lattice of coupled ladders. A vertical line and a dot represent a single ladder and its rung where 
pseudospin Tmn and spin Sm„ (not shown) reside. In the region encircled by the dashed line the Ising couplings between 
pseudospins are indicated. Two pseudospins from (m — l)-th and (m + l)-th ladders are coupled not only by J2 (bold dashed 
line), but also via the dimerization constant e. The 7^„-ordering pattern shown corresponds to the SAF phase, (b): The same 
after mapping on a square lattice. 




FIG. 8: Critical temperature of the PE-SAF phase transition as a function of the Ising coupling g at different values of J, A, e 
from the numerical solution of Egs. 12111 . The dashed line corresponds to the pure IMTF {J = \ = e = Q). Two stars on the 
abscissa show the positions of critical couplings gx = 2.0627 (2.6366) for A = 0.1 (1.0), resp. Large empty circles on the curves 
with £ 7^ indicate the right boundary of the exponential BCS regime 12511 . At large values of g (not shown) all curves T'c(g) 
approach the asymptotic line Tc = g/A. 



